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Abstract 

We investigate divergence- type theories describing the disipative interaction betwen a 
field and a fluid. We look for theories which, under equilibrium conditions, reduce to the 
theory of a Klein-Gordon scalar fleld and a perfect fluid. We show that the requirements of 
causality and positivity of entropy production put non-trivial constraints to the structure 
of the interaction terms. These theories provide a basis for the phenomonological study 
of the reheating period. 

I. INTRODUCTION 

In this paper we investigate divergence-type theories (DTT) describing the dissipative interaction be- 
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I |;ween a field and a fluid [|l],0. We look for theories which, under equilibrium conditions, reduce to the 
Q-i. 

r-| theory of a Klein-Gordon (KG)scalar field and a perfect fluid. We show that the requirements of causality 
. !_1 and positive entropy production put non-trivial constraints on the structure of the interaction terms. 

^ ] The motivation for this work comes from the development of inflationary cosmological models . In 
these models of the Early Universe, most of the energy density is concentrated in a single (fundamental or 
effective) scalar field or "inflaton". There are two well defined moments in the evolution, the roll - down 
period and the reheating period. In the former the inflaton provides an effective cosmological constant 
and supports the superluminal expansion of the Universe in the later the inflaton decays into ordinary 
matter, thus creating entropy and heating up the Universe. 

While inflationary models assume a simple spatial dependence for both the inflaton and the metric of 
the Universe, scalar and gravitational fluctuations around this simple background play an important role 
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1^. During roll - down, scalar and tensor fluctuations are redshifted and frozen by the cosmological expan- 
sion, and become the seeds for future contrasts in the energy density of the Universe. During reheating, 
fluctuations within the horizon interact in a complex way with the background flelds, and are key to such 
questions as what is the duration of the reheating period, at which temperature is equilibrium attained, 
how is energy distributed at the end of reheating, and the spectrum of primordial density fluctuations. 

Overall, quantum fleld theory in curved spaces is an adequate framework to study the roll - down period 
(although of course many questions remain, such as what is the shape of the inflaton potential, and whether 
the flne tuning usually assumed in standard inflationary models is truly unavoidable) [§-|lO|. However, 
reheating is much more complex, as here neither background flelds nor fluctuations may be considered a 
perturbation on the other, and moreover the way the geometry changes along the process has a drastic 
influence on results. For this reason, analysis based on quantum fleld theory on curved spaces have been 
satisfactory only for the early stages of reheating, the so-called preheating where the dominant process is 



particle creation from the background fleld through parametric ampliflcation 1 11-18 



To the best of our knowledge, the most successful strategy to deal with reheating in the strongly 
nonlinear regime has been to describe the inflaton fleld(s) as purely classical [l^. By modelling the 
inflaton as a classical fleld obeying a non linear wave equation it is possible to go beyond perturbation 
theory; moreover, this approach may be justifled on the basis that occupation numbers, in the infrared 
modes of interest, are typically very high. However, a purely classical theory ignores virtual processes 
that may modify the behavior of the quantum theory, as has been shown by recent detailed calculations 
of transport coefficients from quantum fleld theory p0|-|25[|. 

The next step in improving classical models of reheating is to incorporate those quantum aspects as 
phenomenological terms within the classical theory [^. For example, the inflaton and its fluctuations 
may be depicted as a classical dissipative fluid, with constitutive relationships derived in some way from 
quantum fleld theory. 

The simplest and flrst approach to reheating, for example, is to transform the Klein - Gordon equation 
into a telegraphist's equation by adding a dissipative term, where F is estimated from the quantum 
mechanical decay rate of the inflaton ( |27|-p0|). As pointed out by Brandenberger and others this 
approach misses some important components of the reheating process, such as the possibility of preheating, 
namely, enhanced decay through parametric ampliflcation of quantum flelds coupled to the inflaton. A 
subtler criticism is that this approach, or rather, simple covariant generalizations of this approach, are 
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equivalent to writing a first order dissipative theory for the infiaton field (in the classification of Hiscock 
and Lindblom |^), similar to the Eckart theory for ordinary fiuids. However, it is known that first order 
theories have stability and causality problems, and therefore we should expect the same problems will 
occur in these simplified reheating scenarios ( [0). 

The task at hand is then to write down classical dissipative models of reheating, in order to account for 
the basic aspects of the quantum phenomenology, but to do so in the framework of a consistent relativistic 
hydrodynamics, thus building into the model the general principles of relativistic invariance, stability, 
causality and the Second Law of Thermodynamics from scratch. 



In an important precedent to this work, for example, Maartens, Pavon, Zimdahl and others |34-40 



have developed a dissipative model of reheating based on the "truncated" Israel - Stewart framework. 
However, these authors treated the infiaton as a fiuid, thus losing the aspects of the problem associated 
with the coherence of the infiaton field. Besides, they focused on the dissipative effects arising from the 
fact that a mixture of otherwise ideal fluids will generally develop a nonvanishing bulk viscosity. Therefore, 
in their model dissipation happens only if the Universe expands. On the other hand, quantum flelds show 
dissipative behavior even in Minkowski space time. In our work, we shall not only treat the inflaton as a 
coherent fleld, but also focus on dissipation arising from having the inflaton not in equilibrium with the 
ensemble of all other (quantum) flelds. In fact, in this preliminary investigation we shall work on flat space 
time backgrounds, although a fully covariant generalization will be immediate. It is likely that a truly 
realistic account of reheating will require a combination of both this and the Maartens - Pavon - Zimdahl 
approaches. 



Another important difference between this and earlier work is that, instead of the Israel - Stewart |41-43 



framework (or similar frameworks, such as "Extended Thermodynamics" [^) , we shall work within 
the class of "Divergence Type Theories" (DTT) earlier introduced by Geroch |T|]. DTTs are appealing 
because they represent a mathematically consistent, closed system, rather than a truncated expansion in 
deviations from ideal behavior. This makes the discussion of such properties as causality and stability 
most transparent P],^ , for which reason DTTs are a natural language to describe the complex nonlinear 
interactions between the inffaton fleld and the fluid of quantum fluctuations during reheating. 

To summarize, in this paper we explore DTTs describing the interaction between a classical, nonlinear 
scalar fleld and a perfect fluid. Our main goal is to study the constraints on possible interactions which 
follow from the requirements of covariance, causality, stability and the Second Law. For simplicity, we 
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shall work on a flat space - time, and when no loss of generality is involved, we shall assume homogeneous 
configurations. The linearized dynamics of the scalar field around equilibrium is given by the telegraphist's 
equation, but the characteristic dissipative term appears only as a linear and local approximation to a 
more complex term (as it happens in dissipative equations derived from quantum field theory ||^, 



We shall begin our analysis by constructing a DTT whose solutions, under suitable restrictions on 
initial conditions, reduce to solutions of a nonlinear Klein - Gordon (KG) equation. Let us note that while 
it is easy to show that the KG equation is hyperbolic ( [^), to the best of our knowledge this equation 
was never formulated as a DTT. The dynamics of our DTT theory is given in terms of the conservation 
equations for the energy momentum tensor and one vector current ja, and therefore includes one vectorial 
and grange multipliers and ^, respectively. Together with the scalar field 0, these are the 

degrees of freedom of the theory. 

This is in contradistinction with the usual view of and its gradient (p^a as the only degrees of freedom. 
The extra variable is associated to dissipation and will become a true dynamical variable only when 
interaction with other fluids is introduced. The fact that dissipative theories require more degrees of 
freedom than ideal theories is of course rather generic. 

Our next step is to introduce the ensemble of quantum fluctuations of all other matter fields. For 
simplicity, we shall describe this as a single ideal fluid. We shall assume there are no conserved currents 
other than energy - momentum, so this fluid will introduce only one new vectorial degree of freedom j3^^^, 
representing the inverse temperature vector. 

The interacting theory is different from the decoupled theory in two ways. First, the dependence of 
the energy momentum tensors for each component (field and fluid) on the dynamical variables is not the 
same (in a derivation from first principles, this would result from taking the variation of interaction and 



radiative correction terms in the quantum effective action with respect to the metric f^). For simplicity, 
we shall assume that the energy - momentum tensors, but not the current j", are modified in this way. 
Second, only the total energy momentum is conserved, so the divergence of the energy - momentum tensor 
of the field alone, say, is not zero. We may also modify the Klein - Gordon equation (near equilibrium, 
this modification will result in the addition of a F0 dissipative term). As it will turn out, the possible 
modifications in the energy momentum tensors are constrained by causality, while the nonconservation 
terms in the equations of motion are constrained by the Second Law. 

After analyzing the constraints derived from linear stability and causality, we shall conclude this paper 
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by displaying the simplest possible family of theories satisfying all the physical requirements above. We 
shall use these theories for a qualitative discussion of the nonlinear aspects of the approach to equilibrium. 

The rest of the paper is organized as follows. In next section, we show a well defined DTT which 
reduces to the KG equation if suitable boundary conditions are enforced. In Section III we study the 
mixture of field and fluid in the DTT framework. We analyze the equilibrium states and their linear 
stability. In Section IV we write down a simple, acceptable, nonlinear DTT theory, and use it to analyze 
the thermalization process. We summarize the main conclusions in some brief final remarks. 

We gather in the appendix some elementary facts about DTTs and their connection to causality. 



II. DIVERGENCE TYPE THEORY OF THE KLEIN GORDON FIELD 

To investigate properly a relativistic fluid, the context of DTTs proposed by Geroch is particularly 
interesting since causality is easily investigated and/or achieved. We use signature (—,+,+,+), latin 
letters a,6, c, ... indicate spacetime coordinates and Greek letters /i, i/, ... indicates spatial coordinates, 
other conventions follows Misner, Thorne and Wheeler ( [Q). We will also follow the usual (physicist's) 
convention to designate tensorial character by its elements, that is will designate the tensor T = 
T^^Ca ® Cfc, context indicating easily whether one is making a statement about the whole tensor or its 
coordinates. 

Our first goal is to derive a model which contains the classical Klein-Gordon field as particular case. 
Since we want to be able to apply this to the early universe, we allow the classical field to be subject to 
a non-trivial potential (which may even be non-renormalizable pHTOl). The Klein - Gordon theory may 



be described as a conservation law for the "current" 0^a, namely, V°0,a = V'{(f)). The energy - momentum 
tensor is also defined in terms of 0_a as Tab = (p,a<P,b — Qab {4',c4''^/'^ + ^(0))- We therefore postulate that 
our theory is defined by two currents, energy - momentum Tab and a vector current j", with dynamical 
equations 

^"^6 = (1) 

J'^ ■,a = R N (2) 

and the constitutive relation 

Tab ^ ^ajb _ gab f^j^^c ^ y [^A (3) 
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The scalar field cf) is introduced by writing the functional relationship between R and T parametrically 
as -R = V'{(j)), T = V{(f)). This implies no loss of generality. 
Next consider the conservation law 

T'' -fi = 9""' Ub ,c - Jc ,b) f + V if - = (4) 

Then if 

jb ,c - jc ,6 = (5) 

we have = 0'" and we fall back on the usual (classical) Klein-Gordon theory. Note that equation is 
a constraint rather than a dynamical equation. Defining Mat = ja;b ~ jb;a we have the following identity 

V" 

fMab;c = V'Mba + y^f (jbMac - jaMbc) + J^Mbc - j%Mac (6) 

Therefore if (|^) is true initially, it will stay true for all time. In other words, our set of equations (|1|), 
(0) and represent a theory larger than Klein-Gordon, reducing to it if the constraint is enforced 
initially. 

Our next step is to cast this theory within the DTT framework. Since there are two currents, we 



introduce two Lagrange multipliers (we also define /3(c) = ^J~Plc)(^{c)a) as dynamical degrees of 

freedom. ^ is analogous to a chemical potential conjugated to the current j", while /J"^-, plays the role of 
"inverse temperature" and is conjugated to T*^*. Observe the perfect-fiuid form of the energy- momentum 
tensor 

Following the general DTT construction (see Appendix), we introduce the generating function x"" = 
P^^-jP, where p is the pressure and 

^■'^ = /5(c)| (7) 



dp 

r<^b _ ab _ na nb ^ _ ab _ -a -b ^ V / 

Ac) (||) 



leading to the following equation upon comparison with (|^) 

I dp I dp\^ 



P = -ljcf-T (9) 



The first one is a differential equation in the variables ^, /?(c) with solution 



P=^ + ^ (10) 

where A is independent of ^ and /?(c) but can depend on spacetime coordinates. In fact, using the 
second equation we see that 

A = (11) 

Thus 

r = 4rRc) (12) 

P(c) 

Recall that the conservation law for T"^, eq. @) implies — j"0,a = even when the constraint eq. 
(I) is not enforced. Since = —{i/P{c)Y, we get 

^U<t>,a = (13) 



Finally, we easily compute Xc using 



P = - — ^ (14) 



Leading immediately to 



Xe = -^e'lnAc) + ^T/?f,) (15) 

We can regard eq. (0) as a generalization of the canonical momentum vr = 0; in the "rest" frame 
where [3^^^ = (^/3(c),0^ , we get vr = —^/l3{^c)- The problem is that only this ratio has a direct meaning in 
terms of the Klein - Gordon theory alone, that is, the KG equations are invariant under a rescaling of ^ 
and by a common factor. In order to break this indeterminacy, we must look at the larger framework 
where the Klein - Gordon field interacts with other fluids. Then we complete the definition of and /^^^^ 
by demanding that, in equilibrium, (3^^-^ must be identical to the (only) inverse temperature vector of the 
full theory. 

Knowing /J^"^-,, we now regard eq. ([13|) as the definition of the scalar ^. This means that, while for the 
pure KG theory eq. (0) simply follows from energy - momentum conservation, we shall demand it also 
holds unchanged in the interacting theory. This procedure is of course suggested by Landau and Lifschitz' 
treatment of the damped harmonic oscillator in ref. ||51| 

Since xt is a homogeneous function of degree 1 in ^ and (3^^^ the entropy current and entropy creation 
rate vanish in the pure KG theory, as expected for a coherent field. 
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III. DIVERGENCE TYPE THEORY OF INTERACTING FIELDS AND FLUIDS 



Our next goal is to describe the interaction between the KG field and other forms of matter in the 
context of DTTs. We thus introduce a (perfect) fiuid, described by an inverse temperature vector (3^^^ and 
energy - momentum tensor T^^ derivable from a generating functional Xq = 9xq/ df3(q) a (see the Appendix). 
The interacting theory shall be described by the equations 

rpab ja 

-^c \h — 

rpab ja 

-'-q — ~^ 

Ptc)4>,a = -i (16) 

Where j° and T^^ are the current and energy - momentum tensor for the infiaton field, and we have 
added the definition eq. (p^, which, unlike the situation in the noninteracting theory, is now independent 
of the other equations . Let us seek a generating functional of the form 

The total system is generated not only by the sum of each thermodynamic potential but there is a 
third potential to include the interaction between field and fiuid. Each energy - momentum tensor will be 
given by: 



d(3(c) b df3(^c) b 

= + (18) 



Let's define the following variables 



dP{q)b <9/?(g)fe 



= /?f,) - (19) 



Then 



+ 



di3^q)b 2dl3b OB, 



(20) 



Thus 



Tf + T^" = — ^ + -p- + -- (21) 

Note that the only real physical system is the one described by the total energy-momentum tensor 
given above. Any separation in two fluids will be tinged with arbitrariness, a fact that will be clear soon. 
We will ask that the total energy-momentum tensor be symmetric; therefore 

5 will depend in general on scalars as demanded by Lorentz invariance; namely S = E{^,u,v,w, (f)) 
where 

u = -M'' 

w = -PaB" (23) 

Then 

d d d 



dPa du dw 

d d d 

Note that even if the source term P is taken to be null, the energy momentum tensors for fleld and fluid 
do not fall back automatically to their old form. Since the only "true" energy-momentum tensor is the 
total one, it is sometimes helpful and instructive to rewrite the equation of motions of the T^'^ {i — c,q). 
Let us deflne the following 

rjiab rpab i rjiab 

j.ab ^ j.ab _ j.ab ^25) 

Note that 



since dXq/dp^c)b = = dXc/dP(q)b- Also 



rpab^2^ (27) 

oBh 
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The following identities follow straightforwardly 

dE dE 

S« = -2/3»— - S"— (28) 
ou aw 

a'^a ff; a2'^ a2'^ al-^ a2'^ 

^ = -5"'^ + 2/5'^/^"^-^ + 2B-B'-^ + ^(^"B'-^ + B'^(3'^ (30) 
aisfe aw awaw ovow ouov ow^ 

Let 's compute the entropy creation rate. Using 

d d d 

+ 



did Id 



(31) 



we find 

V„5« = ^V„0-S6/^-eA (32) 

where we used the fact that P(^c)^,a + ^ — 0- For simphcity, let us ask that S be independent of 0, so 
that the entropy production reduces to 

VaS" = -Bbl" - (33) 

The second law of thermodynamics imposes that Va-S"" > 0. This means that and A must vanish 
when B'^ and ^ go to zero. In this limit, A, which is a scalar, must take the form A = + Bw {v 
being of higher order); Lorentz invariance demands l'^ — — {C^ -\- Dw) ff' — EB^. Therefore Va^*" — 
_(^B + C)w^- Dw'^ - Ev, and we must have A, D < 0, E < and AD > {B + Cf /4. Equivalently, 

we may parametrize A = M««, Bp" = 2(1 - k)M^'', C(3^ = 2kM^^ and L'/3"/3'' + Eg"^ = M'^^ whereby 

A = M^^{ + 2Bail - k)M^" (34) 
= 2k^M^'' + BaM""^ 



A. First order analysis away from equilibrium in a simplified, homogeneous model 



Let 's turn now to the equation of motion to analyze small deviations from equilibrium. The requirement 
that the equilibrium state must be isotropic in some "rest" frame implies that, in equilibrium, the two 
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vector and must be parallel. Therefore we can write jS'^^-^ = jScU"" and ^ = with a common 
unit vector m". We will also restrict ourselves to the homogeneous case. Observe that T^l = T^\^ = 
identically for /i 7^ 0, so there are only four nontrivial equations, including eq. (p!3|). 

In equilibrium, we must have ^ = A = I^ = R = 0. Defining = at the equilibrium point, we conclude 
i^eq , Peq, Beq , (peg) = (0,/?o,0,0). We cau uow analyze the first order deviations away from equilibrium; 
that is 

e = 5e 
(f) = 6(f) 

I3 = l3eq + 513 

B = 6B (35) 

We can expand 

00 

E = Y.Fn{P,OB-^' (36) 

n=0 

H is constrained to this form by demanding that at equilibrium, where f3c = Pq, we fall back to our 
initial T"'' and j". A and J* are given by eq. (^4|), where all the matrix element are evaluated at the 
equilibrium value ^eq = = 0, = (3l^ and 0o = 0. Moreover we will define w? = V" (0), so to first 
order 

R = m'^S(j) (37) 

In our case we have 

j% = ^ {SO, = Ml^S^ + 2(1 - ^)Ml'SB + m'S<P 



d/3dB^ 

= p'q m,t + 



Po (50) ,i = -5^ (38) 
We can revert to our old coordinates SPc and 6Pg 

^ (50 t = M^^S^ + 2(1 - «:)M«°(5/3e + 2(1 - K)M^o%(3q + mH(P 
Po 
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dpdB^ 



d/3dB^ 



The advantage to do this is that we can add the middle two to obtain 



That is 







S(3q — constant 



We are thus left with 



or, in matrix form 



^ m-' 





\ 













V° 




J 





M^^6^ + 2(1 - k)M^%P^ + 



^ M^^ 2(1 - 
2«;M«° MOO 

V-1 



2 

m 


\ 














/ 





We will consider solutions of the type exp {—jt) . There is a solution if 



(-7/9o-^ - M««) -2(1-«:)M, 



CO 



—m 







-MOO) 



-7/3o 



= 



To make notation less clumsy, let us write 



d(3dB'^ 



Mf = 



m; 



00 



si^ A 
X 

£3 = Dl3l + E 
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The notation is chosen to emphasize exphcitly which quantities are small and which aren't . Namely, 
the Ei are (very) small but the X is not. Therefore 



(-7/3o-'-£i) -2(l-«)£2 -m^ 
-2Ke2 (7^ -£3) =0 

1 -7/3o 

Expanding with the second column we have 

The solutions to order zero are 7 = ±im and 7 = 0. To order one, the equation reduces to 

XPo-f^Si + X-f (72 + m^) - £3 (7^ + m^) = 



(46) 



(47) 



(48) 



Writing first 7 = ±im + S we obtain 



(49) 



and in the case — S we have 



(50) 



That is 



7 = ±im 
£3 

^^X 



(51) 



Notice that £1 < and £3 < . Since Po > then the first two solutions corresponding to oscillatory 
modes get a damping part. The other solution will also be damped if 



X = < 



(52) 



dpdB^ 

In this case the system will be causal if the matrix M^'^^v is negative-definite (where f = f " > and 
vq — —v^, since it is the temporal component of a temporal 4- vector oriented toward the future) and where 

^ vpo^ ^ 








-v(3o 
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Obviously, this matrix is negative-definite if 



/?o>0 



< 



d(3dB^ 

We tlius verify that the causahty condition imphes the stabihty of the equihbrium solution. This is in 



agreement with the findings of Hiscock and Lindblom in the context of the Israel - Stewart formalism 
Going one order further we find 



7 = ±i < m — 



£3 



/3o_ 
2m 



1^^2 4k (1 - k) 2 
-fJo£i ^ £2 



(53) 
(54) 



That is the correction to the third solution is at least third order and the second order terms modify 
the frequency of the oscillation in the case of the other two solutions. 
If we look back to the linearized equation of motion for 6(j), we obtain 



- (50) = -M««/3o {6<P) , + 2(1 - K)M',''Sf3, + m 



(55) 



In the simplest case Mq^ = 0, this is the telegraphist equation with a damping term Tcp, where 
r = Mq^ Pq. It is important to realize, however, that this identification holds only to linear order away 
from equilibrium. In general, F will not be a constant, but a function of the dynamical variables f]q and 
(3c, and therefore it will depend on the history of the system. 

We may observe that when Mq" = 0, we get Pc,t = to first order. 



IV. A SIMPLE NONLINEAR MODEL OF FIELD - FLUID INTERACTION 

In this section we shall investigate the simplest DTTs of field - fluid interaction satisfying the re- 
quirements of causality, stability and the Second Law discussed in the previous Section. To make things 
simplest, we will restrict the form of S in the following manner 

S = F{u)v + G{u)w'^ (56) 

and we will investigate what restrictions one can expect from causality on F and G. By making S 
independent of ^, we make sure that the current j'* preserves its form in the interacting theory. We also 
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assume S to be cj)— independent, and require it to be quadratic on the difference variable 5", in such a 
way that corrections to the energy - momentum tensors will vanish in the equilibrium state. Observe that 



= -2(3'' (f'v + G'w^) - 2B''Gw 
4/3"/3^ [f"v + G"w^) - 2g''^ [f'v + G' 



+4 {^''B^ + G'w + 2B''B''G 



Ap^'P^G'w - 2g'''Gw + A^^'BT' + 2B''(3'G 



^ab 



(57) 



At equilibrium 



Q2^a 



(58) 



and 



OB.dB, 
Causality demands that 



(59) 



(eg) " g 



0^ 



V 


















(60) 



J 

should be negative-definite for any future-oriented, timelike vector v"" . Since Xq represents a perfect 



fluid, we obtain causality under the usual conditions [|1| . We only need to verify that g^^^Va is negative- 
definite by itself. Specifically, we want to know if 



■Ja = AP^eq)Va (-/9fe,)/?(e,)G" + g'^Ff) + 2 {v' ^l^^^ + V^^l^^) G 



(61) 



OB^dB, 

is negative-definite, that is if Li,Lc-^^^Va < for any La- To achieve this, let's decompose Lb into its 
longitudinal and transverse part relative to P{eq)b '■ 



U = l 



We easily work out the following 



u 



+ Rb 



(62) 
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(63) 



We can already extract some information. Taking / = we have 



(64) 



which imply that F' > since P^^g^Va < and RbR'' > . Also, Taking i?;, = we obtain 



U 



(65) 



this time implying that i—G' 



(F'+G) 



> 0. Let us now write 



(66) 



and decompose the spatial vector Ri, into the part which is longitudinal and transversal to a;*: 



R^ = r]uj^ + 



(67) 



Therefore R^v'^ = rjUhU^ and, replacing in (|6lD , we obtain upon division by /^^ggj^a < 0: 

which should be now positive-definite in order to have (|6TD negative-definite. Since Ut^F' > it suffice 
then to ask 



-G 



, {F' + G)\ u;,u;' , ^2^,^,.^, 



u 



-Gr]l + r]'ubu''F' > 



(69) 



This is a quadratic form in / and 77 and we are thus demanding that the matrix 

/ / ^/ {F'+G) \ u^.u," ^\ 



(70) 



should be positive-definite. We already know that this means that the diagonal element are positive, a 
fact that we already deduced. Now since the matrix is real and symmetric, we know that its eigenvalues 
are real. Since the diagonal elements are positive, then it is sufficient to prove that the determinant is 
positive. The determinant is given by 



u 



4 I Pl,gfa 



> 



(71) 
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This condition can be simplified by considering that VaV'^ = ^^P(eq)aP(eq) + ^a^^"" < that is, using 

LUaUj'' < \\ (72) 

Also P'^eq)'^a = —\u . Therefore 

^a^° _ UgU" 1 
2 ~ \2,,4 ^ 



(eg) 



Va 



Then if 



f_G'_(£:±^V'>^i (74) 
then condition (fflD is fulfilled. That is, the matrix (??) is positive definite if, 

F' > (75) 

{G'^^£±^y, (76) 

and 

\ u ) Au 

To see the meaning of these conditions, consider the case where F and G follow power laws. Dimensional 
analysis in natural units {c = h = ks = 1) indicates that 

[H] = (78) 

Thus, writing F = Fqu"'^^ and G = Gqu" implies that a = —3 , that is 

S = Fqu-^v + Gou-^w^ (79) 

The three conditions (??),(??) and (??) read 

Fo<0 

Go + Fo > ^ Go > -Fq > 

AFoiGo + Fo)<-^Gl (80) 
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Substituting Fq/Gq= —x, we have the following inequalities for x 



- 4x + - < 
4 

implying that 1/2 — -\/3/4 < x < 1/2 + •\/3/4. Observe that both hmit cases x 
excluded. 

We thus obtain the following final form for S 



(81) 



and X — 1 are 



u 



X 



u 



(82) 



(f]o 



— Go < 0) which can be rewritten in a covariant manner as (writing ji = v^^3a^) 



5 = ^B^B' 



-xgab 



(83) 



A. Nonlinear approach to equilibrium in the homogeneous case 

We shall conclude this paper by studying the nonlinear evolution of the simple model just described. 
By simplicity, we shall work in fiat space time and assume a homogeneous model, implying in particular 
that and B'^ are colinear, = and B"' = Bu"'. We choose a frame at rest with the fluid where 

The equations of motion are 



jl = F'(0) + A 



(84) 



(85) 



',0, 



(86) 



Lets introduce the two energy-momentum tensor of both (perfect) fluids without interaction 



(87) 



-V{<P) 



and 
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The interaction between the two is given by (/3 = ^{Pc + Pg) and B — 

rpab _ 



and 



rpab 2 



We use as model for our interaction term 

Under the simphfying hypothesis stated above we have 



r. 



2 /32 



00 



+ V{<t>) 



and 



as" _ B_ 



It is convenient to introduce a new variable s defined as 



implying 



Pc + Pg 



and to express the equations in terms of the canonical momentum 



Pc 



The energy-momentum tensors thus read 



TT 



a + 2or- 



1 + s/2y\ P 
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71 



T^"^:- + V (0) - 



CT + ler 



{l + s/2f 



1 



We are left with four O.D.E. 



-7r,t = V'{<P) + A 

(r°°) ^ = 2/° 



where 



and 



7° = -SM°° + 2<M^° 
The dimensions of the M^^ terms are 





^2 








rp4 





This vahdates the following modelling for the M 



AB 



where r is some function of the dynamical variables with dimensions T ^ and Kq^Lq > 0. Moreover, 
positive entropy production imposes the following relationship: 
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KoLo > Ml 

We will also take the q -fluid as describing radiation so that p = aT^ with the Stefan - Boltzmann constant 
cr > 0. The equation (|8^ becomes 

-vr, = V'{<P) + ^ (s + 1) /3,7r - 2 (1 - k) sf3,^ 

Observe that 



And define 



y + V ((/)) 



TT [n^t + V 
Ko 



TT 



S + 1) PgTT + 25 (1 - K 



Mn 



/° = Lo^/3,-/€(s + l)§7rMo 



Mo 



ci = a + ler 



(l + ./2)^ 



Then 



and 



C2 = a + 20r 



(1 + ^/2)' 



7r,t = -V (0) - K 



4 ^ 40r sil - s) 



U -2/0 



The final system of equations is 



/3,ni+g/2/ 

8rc3 



{2/% + ir7r (C1 + C2)} 
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TT,, = -V'{cf>) - K 



TT 



where 



V (0) = ^m^' + 



Finally 



4 



is both positive definite and conserved. 

These equations generally describe the approach to equilibrium. This is most clearly seen in the limit 
where s remains small. In this limit, we have Ci ~ C2 ~ cr, C3 ~ —2a < 0, /i (s) ~ 2 and K ~ K^j^qi^jT^ . In 
particular, the equation for the field describes damped oscillations, but the damping "constant" is truly 
a dynamical variable, thus opening a mechanism to include memory effects in the dynamics and interesting 
behavior; note by exemple that K can change sign if the coupling term Mq 7^ 0. 

More generally, by exploring parameter space we find a variety of behaviors. In particular, the approach 
to equilibrium may be either over or underdamped; this second case seems to be relevant to the description 
of preheating episodes. 



V. FINAL REMARKS 

The main conclusion of this work is that the conditions of causality, stability and a proper thermody- 
namic behavior put concrete limits on possible phenomenological models of the reheating period. We have 
shown concrete examples of divergence type theories which satisfy these requirements. Unlike earlier work, 
we have described the inflaton as a field, rather than disregarding its coherence by describing it simply as 
another fluid. This has required an extension both of the usual Klein - Gordon and DTT frameworks. We 
have also shown how this field can be consistently coupled to a fiuid. 
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The equations of motion we have derived in this last Section show also that it is possible to explore the 
bulk of possible dynamical behaviors already with models with a minimal set of undetermined parameters 
(in our case, these were F, Kq, Lq, Mq and the functional form of r). These parameters may be estimated 
by fitting the predictions of the model to microscopic calculations in controlled limiting cases, much in 
the same way as viscosity coefficients in field theory are computed by analyzing the damping of extra long 
wavelength fiuctuations p0|-[25[]. We may then obtain reliable phenomenological models to use as a tool to 
explore the full nonlinear physics of reheating, with an enormous gain in simplicity as compared to a full 
attack from a first principles perspective. 



VI. APPENDIX: DIVERGENCE TYPE THEORIES 

Following Geroch divergence type theories are usually described in terms of some tensorial quantities 
that obey conservation equations 

rj^ab ^ Q 

A^-a = 

j^abc jbc 

;a 

This is a simple and slight generalization of relativistic fiuid theories proposed initially by Liu, Muller 
and Ruggeri [^]. In this setting, T"''' is the energy-momentum tensor and N"- is the particle current. Their 
corresponding equation simply expresses conservation of energy, momentum and mass. The third equation 
will describe the dissipative part. The energy- momentum tensor is symmetric and A°''"^ = A"''*\ A"-^ 6 = 
and = 0. The entropy current is enlarged to read 

The ^, ^"'^ are the dynamical degrees of freedom. The following relations hold 



T 



ab 



A 



abc 



dihc 
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Symmetry of the energy-momentum tensor implies that 

That is all the fundamental tensors of the theory can be obtained from the generating functional x- 
The entropy production is given by 

Positive entropy production is ensured by demanding that I^^ = M^'^'^^^'^'^^^de, where M is negative 
definite. 

Ideal fluids are an important if somewhat trivial example. To obtain ideal hydrodynamics within the 
DTT framework, consider a generating functional Xp = Xp {^^ f^) where fi = \/—^^- It is a simple matter 
to obtain 

e dxp 



/i djj, 



^ ji dji jj? n djji djj? 

A simple comparison with the perfect fluid form of the energy-momentum tensor T"'^ = g"'^p + 
y^ayb ijjiplies the following identification 

P = --^ (88) 



P 



Note that the conserved current can be quite generally written as 

A less trivial but important example both historically and conceptually is the Eckart theory which can 
be obtained from 

XE = Xp+ ]^CabU''v!' 

Performing a Legendre transform to the new variables one obtains a system of first order 

differential equations of the form 
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where stand for the entire collection of variables (^,^a,^afe) and similarly = (^0,0,/"''^ represent 
the dissipative source; the index A thus covers 14 dimensions in our example. This first order system of 
differential equations is symmetric since 

Note that we have a system of the form 

A'v^i + Bv = 

where z is a space-time index, the and B are k x k matrices and v is a k -vector. Now this (first 
order) system is hyperbolic if all its eigenvalues are real; each of these eigenvalues represent the velocity 
of propagation of some small disturbance in space. These in turn propagate along hypersurfaces called 
characteristics whose existence is insured by the existence of k real eigenvalues []53| , ^ . If the matrices 
A^ and B are symmetric then it suffices that some combination A^Vi be definite (negative-definite given 
our choice of the signature for the metric) to insure that all the eigenvalues are real (but some could be 
degenerate). An usual case happens when this combination reduces to A^ , the vector v being the time- 
like vector (l, if ^ . In a relativistic theory one would expect hyperbolicity to be invariant under (proper) 
Lorentz transformations; in this case we say the system is causal. In our context, one would thus say that 
the system is hyperbolic if 

Q2^a 

is negative-definite for some temporal vector Va and the theory will be causal if this stay true for any 
temporal vector Va- 
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